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Abstract 

Generalized parton distributions provide a unified parameterization of hadron structure and allow one 
to combine information from many different observables. Lattice QCD calculations already provide 
important input to determine these distributions and hold the promise to become even more important 
in the future. To this end, a reliable extrapolation of lattice calculations to the physical quark and 
pion masses is needed. We present an analysis for the moments of generalized parton distributions of 
the pion in one-loop order of chiral perturbation theory. 



1 Introduction 



Generalized parton distributions (GPDs) provide general parameterizations for non-forward hadronic 
matrix elements, see e.g. the original work PH2HHJ and the recent reviews lUISHS]- They contain the 
usual parton distribution functions and elastic form factors, and they are closely related with light- 
cone wave functions and distribution amplitudes. GPDs can be analyzed using standard operator 
product expansion techniques |HE|, and factorization theorems have been proven for many different 
processes. In addition, evolution equations have been derived to next-to-leading order accuracy, 
and generally speaking the GPD formalism has reached a similar level of stringency as the QCD 
description of inclusive deep inelastic scattering. Finally, it has become clear that GPDs contain 
information which cannot be deduced directly from any experiment, most notably the orbital angular 
momentum carried by partons [2] and the transverse structure of hadrons in the impact parameter 
plane |Hj. On the experimental side, data which is especially suited to constrain the form of GPDs 
has been obtained at DESY |Hj and Jefferson Lab |lUj . 

An important problem hampering the practical realization of the physics potential of GPDs is 
that, typically, measurements are only sensitive to convolutions of GPDs, such that the assumed 
functional form enters crucially into the extraction of the distributions themselves. The situation is 
further complicated by potentially large higher-twist contributions in some processes. Therefore, the 
possibility to obtain independent information from lattice QCD, which allows one to directly calculate 
x-moments of GPDs, is especially interesting. Pioneering lattice calculations for nucleon GPDs have 
already been performed some time ago These calculations are being improved and their scope is 
being systematically widened |12j . also in the direction of analyzing other hadrons, in particular the 
pion. Calculations for the pion are being performed by the QCDSF collaboration, using dynamical 
Np = 2 improved Wilson fermions, but have to be made at unphysically large quark masses to save on 
computing time. The reliable extrapolation to the chiral limit, i.e. to physical quark masses, requires 
precise knowledge on the functional form of the quark mass dependence of moments of GPDs. Close 
to the chiral limit, this knowledge is provided by chiral perturbation theory (ChPT), and in recent 
years the combination of ChPT and lattice QCD has in fact proven to be very powerful in obtaining 
precise lattice predictions. Obviously the ChPT analysis should be performed in parallel to the lattice 
calculations, and this is the objective of our paper. 

In this Letter we present an analysis for the moments of pion GPDs in one-loop order of ChPT. 
The paper is organized as follows. In Section [2] we recall the definition of the pion GPDs and the 
expression of their moments in terms of the matrix elements of local twist-two operators. In Sections^ 
and 0] we construct the relevant twist- two operators in the first two orders of ChPT, and in Section [5] 
we calculate the one-loop contributions to their matrix elements in the pion. We summarize our 
findings in Section HO 

2 Definitions and properties of pion GPDs 

The pion GPDs can be introduced as the pion matrix elements of nonlocal operators. As is done in 
lattice QCD calculations, we assume isospin symmetry to be exact, neglecting the difference between 
u- and d-quark masses. We shall study the isoscalar and isovector quark GPDs, which are respectively 
defined as 
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where q denotes the doublet of u- and d-quark fields and r a the Pauli matrices in isospin space, u is 
a lightlike auxiliary vector, and we use the standard notations for the kinematical variables 

d P + P' a / + A 2 c Au (?\ 
P = , A = p — p, t = A , f = — — . 3 

In terms of the distributions for individual quark flavors defined in [£] we have H I=0 = H^ + H% and 
H I=1 = H% — H!~. As usual, Wilson lines between the quark fields are to be inserted in (JTJ and (j2J) 
if one is not working in the light-cone gauge u^A^ = 0. 

Time reversal and charge conjugation symmetry result in the following properties for the pion 
GPDs: 

H\x,i,t) = H^x-^t), (4) 
H\x,^t) = {-if-'H^-x,^). (5) 

The distributions can be analytically continued to positive t and are real- valued for t < 4m^, as 
discussed in ^3 ■ 

The x-moments of the GPDs are related to the matrix elements of the local isoscalar and isovector 
operators 

1 «-» «-» 1 

where = ^(-D^ — D^) = ^(8^ — 8^) — igA^ and all fields are to be taken at the same space-time 
coordinate. S denotes symmetrization and subtraction of trace terms for all uncontracted Lorentz 
indices. Defining form factors 

n 

Wl^-^KW) = §ab E ^WSV-.A^...?^, (7) 

fc=0,2,... 
n 

(*\p')\O c ^ 2 ..., n K{p)) = ie abc A I = k 1 (t)SA, 1 ...A, k P, k+1 ...P, n , (8) 

k=0,2,... 

where the operators are to be taken at zero space-time position and where time reversal invariance 
limits k to even values, one finds 

-1 n 

/ dxx"- 1 H I (x,£,t)= E <*(*)• (9) 

1 fc=0,2,... 

Note that because of the symmetry relation (J5J the moments with odd n in Q vanish for 1 = and 
those with even n vanish for 1 = 1. The form factors A^ k (t) are real- valued in the region t < 4m 2 . 
Information on the lowest moments 

J dxH I=1 (x,£,t) = A{f(t), (10) 

j\xxH I=0 (x,U) = Aif^ + AfAifit) (11) 

will soon be accessible from the lattice for pion masses down to values m n ~ 500 MeV |15j . Note that 
the first moment of H I=1 is related to the pion electromagnetic form factor as Aj^(t) = 2F£ m (t) 
and the second moment of H I=0 with the form factors in the parameterization ^H] 



(Ap')\T®\K a (p)) = 6 ab 



2P„P V 2 (i) + \ UluA 2 - A,Au) 9i(t) 
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where T® u = \q (7^ iD v + 7^ i-D^) g is the quark part of the energy-momentum tensor. Comparison 

with Q gives 2 (*) = 4£o°(t) and °i( t ) = - 4j4 2^°(0- 

An analogous discussion can be given for the gluon GPD, which is defined by 

= ±J ^ e »W (ttVJI^^C-^^/C^Au) |tt») (13) 

in the convention of 0- It has the same symmetry properties as xH I=0 (x,^,t), and its moments 
/ dxx n ~ 2 H^.(x,S,,t) are analogs of the moments J dx x n ~ 1 H I=0 (x,S,,t), with which they mix under 
renormalization scale evolution. For ease of notation we will in the following concentrate on the quark 
sector, bearing in mind that our results for the isoscalar form factors A^r^{t) readily generalize to 
their gluonic counterparts. 



3 Tensor operators in ChPT 

We use the standard 0(p n ) power-counting of ChPT, where p denotes a generic pion energy or 
momentum. The Goldstone boson fields are collected in the matrix-valued field U, and the leading- 
order Lagrangian, which is of order 0(p 2 ), reads 

4 2 2 = ±F 2 tt(p ft UD' t rt + x i U + xU 1 ) (14) 

for two light flavors, with the covariant derivative and the external field tensor x defined as 

DfjJJ = d^U - ir^U + iU , (15) 
X = 2B (s + ip) . (16) 

Here s and p denote external scalar and pseudoscalar fields, which count as quantities of order 0(p 2 ). 
Furthermore, and L are external right- and left-handed vector fields with intrinsic chiral power 
0(p). The two leading-order parameters appearing in (|14|) are the pion decay constant F (normalized 
to F ~ 92 MeV) and the two-flavor chiral condensate Bq = — (0\qq\0)/F 2 , both evaluated in the chiral 
limit Throughout this work we use the non- linear representation for the pion fields 7r a , 

U = exp{iT a TT a /F} . (17) 

Left and right chiral rotations of the nonlinear pion field 

U^V R UVl (18) 

induce the following transformation of the external fields: 

r, - V^V^ + iV^Vl x ^ VrX vI (19) 

One method to calculate the matrix elements of the operators O s and „ from © 

in ChPT would be to introduce corresponding tensor currents as external sources in the Lagrangian, 
in generalization of the vector currents and L in (|14|). For tensors of higher rank, this becomes 
however increasingly awkward. Instead, we will directly match the operators in onto operators 
with pion degrees of freedom. The expression of the isovector vector current is for instance well-known 
to be plj 

Ol = \qi^r a q ~ V£ = --F 2 tr [r a (rfdJJ + Udjfl) } (20) 
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to leading order in the chiral expansion. Despite the progress in deriving chiral dynamics directly 
from QCD (see e.g. ^Sj) this matching problem is far from being solved in general. The so-called 
low-energy constants appearing in ChPT can in principle be determined by fits to either experimental 
or lattice data, but only a limited number of them are known in practice. What we can however do 
is to construct for a given quark or gluon operator in QCD all operators in ChPT with the same 
symmetry. Their matching coefficients typically are linear combinations of low-energy constants and 
can be fitted directly to lattice data at different unphysical pion masses, provided they are sufficiently 
close to the chiral limit. In Section [5] we will determine the functional dependence on m\ of the form 
factors k (t) required for such a fit at one-loop level in ChPT. 

For the discussion of chiral symmetry it is useful to consider the left- and right-handed tensor 
operators 

(0§W../* = Sq ]1 ^J—^iD^...iD lin q u (21) 

d 1 + / Yk *-* <-> 

(Og) w .„^ = Sq j ^—^-iD w ...iD lln q i , ( 22 ) 
which take values in two-dimensional isospin space. The operators in © are then given by 

<^...Mn = \^ + °t..,n } . = ^tr [r a (0^ n + 0* } • (23) 

Note that for ease of notation we shall use the same symbol for the quark operators in QCD and for 
the operators representing them in ChPT. The left- and right-handed operators behave like 

Ojx..^. - V L L ^VI 0*^ n - V r O*^V r (24) 

under chiral rotations and are related as Of!' „ = PO^ ,, P~ l by the parity transformation P. 
To construct a linearly independent set of operators for the expansion of 0% 1---fin , let us inspect the 
building blocks at hand. For the purpose of calculating matrix elements between pion states we 
can set the external fields p, and L to zero, and keep the scalar field s as a nonzero constant 
field implementing the explicit breaking of chiral symmetry by the quark masses. Since we assume 
exact isospin symmetry, this simply gives x = m2 lj where m is the bare pion mass in ChPT. For 
constructing the operator this leaves us with the objects U, x an d d^, with the additional condition 
dfjX = 0- ^ is convenient to introduce the following quantities 

L^tfd^U, x± = X ] U±U^x, (25) 

which transform covariantly under left rotations and are invariant under right rotations. The product 
of some number of and x± with any number of additional derivatives acting on them transforms 
covariantly under left rotations. Since d^x± = \ [x±>-k/u] + \\x%^ji\ f° r d^X = 0) one can eliminate 
all derivatives acting on x±- Similarly, taking into account the identities 

d^L u - d u L^ = [L^,L U ], dpdPL v = $,0%, + WL V ,L^\ + (26) 

and the equation of motion for the field U, 

d*L lt = -~(2 J £-ltr X t), (27) 
one can trade the tensors d Ul . . . d Un _ 1 L Un for the symmetric and traceless ones, 

L Vl ...v ri = Sd Ul ... d Vn _ x L Vn . (28) 
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A linearly independent set of operators transforming like C^..- under chiral rotations can hence 
be constructed from products of the fields x± an d the tensors L ul ___ Uk . In the ChPT expansion of 
^ n they can be multiplied by chirally invariant operators, which can in turn be constructed from 
the traces of products of x± and L Vl ^ Mk } In each term of the expansion, some Lorentz indices of the 
tensors L Vl __ Vk are equal to fii, . . . ,fj, n and all others are contracted. Since <9 M and are of order 
0{p), each contraction of the indices, as well as the insertion of a field x± increases the chiral counting 
by a factor 0(p 2 ). 

At leading order of ChPT the operator „ n can hence be represented as the sum of products 
of L tensors and their traces, with no contracted Lorentz indices. It is convenient to contract free 
indices with the light-line vector u, introducing the notation 

O h n {u) = OjL^u" 1 • • • L *( u ) = Wx-^ 1 ■ ■ • (29) 

The operator 0^(u) is then expanded on terms of the form 

L kl (u) . . . L kji (u) tr { L kji+1 («)... L kj2 (u) } tr { L kj2+1 («)... L kja (u) } . . . (30) 
and explicitly reads 

°n( U ) = c nL n {u) + ^2 c n,k L k {u) L n _ k (u) + . . . , (31) 
k=l 

where we used that tr = and have omitted terms containing three or more L tensors. Such terms 
do in fact not contribute to the two-pion matrix elements in leading order and next-to-leading order 
of ChPT. To see this, notice that the expansion of into pion fields starts as 

^ = ir a (d^ a /F - e abc 7r b d^ c /F 2 + |vr fe {v a d^ b - ir b d^ a )/F s + 0(vr 3 ) d^) . (32) 

Terms with three or more L tensors hence do not contribute to two-pion matrix elements at Born 
level. At one loop, they contribute only through the tadpole diagram of Fig. where three or four 
pion fields from the inserted operator then carry a derivative contracted with the light-like vector u. 
The numerator in the corresponding Feynman integral thus has the form (uk) 1 , where i > 1 and k 
is the loop momentum. For odd i the integral vanishes because the integrand is antisymmetric in k, 
whereas for even i the integral is proportional to (u 2 ) 4 / 2 by Lorentz invariance, and hence zero as 
well. 

Repeating our discussion for the right-handed operators OS _ „ , with 

R lt = Ud lt U^, x± = xU ] ±U X ] (33) 

and definitions analogous to (|2*5|) and (J2HJ), the tensor operator O n {u) = 0^{u) + 0^{u) reads at 
order O(p ) in ChPT 

n-l 

O n (u) = c n \ L n (u) + R n (ufj + ^2 °n,k \ L k (u) L n _ k (u) + R k {u) R n _ k (u)j + . . . , (34) 

k=l 

where the dots denote terms which do not contribute to two-pion matrix elements to 0(p 2 ) accuracy. 
The coefficients of the L and R terms must be equal because of parity invariance. Using the hermiticity 

1 Equivalently, one may take the traces of products of \± an d Rv x ...v k introduced below. 
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and time reversal properties of O n (u) following from (|21j) . (j22j) and the corresponding properties of 
L n (u) and R n (u), one finds 

(Cn)* = -Cn, (C n )* = (-l) n C n , (35) 

(c n ,k)* = Cn,n-fc , (c n ,fc)* = (-1)™ C„,fc . (36) 

Neglecting terms which do not contribute to two-pion matrix elements one derives for the isoscalar 
and isovector operators at order O(p ) 

n/2 

O s n (u) = j2a n>k (L a k (u)L a n _ k (u)+Rt(u)R a n _ k (u)), (n even) (37) 

k=l 



O a n (u) = ib nfl (L a n (u) + R a n (u 



(n-l)/2 



+ e abc Yl b ntk (L b k (u)L c n _ k (u)+R b k (u)R c n _ k ( U )), (n odd) (38) 

k=l 

where L k (u) is defined by L k (u) = r a L k (u) and R k (u) in analogy. The expansion coefficients a n ^ k 
and b n)k are linear combinations of the coefficients in (|34j) . 

a n ,k = c n ,k + c n>n _ fc for k > 1, (39) 

b n ,k = «(c n ,fc - c n) n- k ) for fc > 1, 6„ i0 = -ic n . (40) 

Because of the symmetry relations (|35j) . (|36JI one has a nj /~ = for odd n and b n ^ k = for even n, 
and the nonzero coefficients a Hik and 6 n & are real. To obtain simple expressions for the form factors 
k (t) in ©, we rearrange derivatives and obtain for the isosinglet operator 

n-2 _ 

O s n (u)=F 2 E <k(iud) k ( K L a {u){2iud) n - k - 2 L a {u) + R a {u)(2iud) n - k - 2 R a {u)) (41) 

k=0,2,... 

with d = r;(d —d). Here we used L k (u) = (ud) k ~ 1 Lf(u) and the abbreviation L a (u) = Lf(u) = u^L® 
and their analogs for the right-handed fields. The coefficients a n ^ k are linear combinations of the a n>k . 
Calculating pion matrix elements one gets for the isoscalar form factors A^jPft) at order O(p ) 



^(0) = 2 n ~ k 



0"n,k-2 — a n,k 



(n and k even, k < n) (42) 



where we set a n ^ n = a n -2 = 0. To ease the notation we have omitted the isospin label, implying here 
and in the following / = if the first index n of a form factor is even, and I = 1 if it is odd. The 
inverse of (|42|) reads 

k n 

~a n>k = - E 2 m - n 4°L = E 2 m " n 4°L- (n, m and k even, k < n - 2) (43) 

m=0 m=k+2 

Prom (J42[) we also see that at order O(p ) the sum of the isoscalar form factors with given n vanishes, 

n 

E 2fc <l = > (™ even ) ( 44 ) 

fc=0,2,... 
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which according to © gives the soft pion theorem H I=0 (x, £ -► 1, 0) = first derived in [TI]. Recalling 
our discussion at the end of Section |5J we emphasize at this point that the ChPT equivalent of the 
twist-two gluon operators has the same form as Q41JI . of course with different matching coefficients. 
Analogs of (|42j) to (|44[) therefore hold for the appropriate moments of the gluon GPD H^{x,^,t). 
Rewriting (|38|) in a form similar to (|41|) gives for the isovector operator 



O a n {u) = b n>n ^(iud) n - x V a {u) 

n— 3 

+ iF 2 e abc Kk (iud) k (L b {u) (2iud) n - k - 2 L c (u) + R b {u) (2iu3 ) n - k - 2 R c (u)) , (45) 

fe=0,2 

where the b n ^ are linear combinations of the b n ^ in p8|) . Here 

V a (u) = u»V* = - l -F 2 (L a {u) + R a {u)) (46) 

is the isovector vector current already given in (|2Uj) . which also implies b\ o = 1. For the isovector 
form factors A^{t) one derives at O(p ) 



A^l = 2 n ~ k 



k 

^n^k — b n .k-2 



b n , k = 2m ""4 0) m> (n odd, k even, fc < n — 1) (47) 

m=0,2,... 

where we set 6 n ._2 = 0. Note that the isovector analog of (|H)) . 

n-l 

2 2 fe 4°i = 2 n 6 n , n _ 1 , (nodd) (48) 

fe=0,2,... 

does not have a vanishing right-hand side. The coefficient 6n,n-i is however related to another 
observable, namely the one-pion matrix element of the isovector axial twist-two operator 

Oln{u) = \tr [r a (0*(u) - O^u)) } = \qu^ (iuD)^ 1 r a q , (49) 

which to leading order of ChPT has the form 

Ol n (u) = bn^iudT^A^u) + ..., (50) 

where A a (u) = u^A 1 ^ = —^iF 2 [R a (u) — L a (u)] is the isovector axial current and the dots denote terms 
which do not contribute to the one-pion matrix element to 0{p 2 ) accuracy for reasons analogous to 
those we discussed after (J22J). One readily finds (vr 6 (p)| 0% n {u) |0) = -iS ab (pu) n F 

bn,n—i a t order 

0(p ). The pion distribution amplitude is defined as 

i/g e" A(P " )/2 (* h (p)\ ?(-5A«) u» W r a q{\\u) |0> = -iS ab F n <j>„{x) , (51) 

where ^(1 + x) is the momentum fraction of the quark in the pion and i J —1 cte 4> n (x) = 1. Here 
is the physical value of the pion decay constant, given in one-loop ChPT as ^7] 

^ = F (l - log ^ + ^ (JO*) + eV)) , (52) 
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where \i denotes the renormalization scale and is a renormalized low-energy constant. Taking 
x-moments of the distribution amplitude, one obtains the matrix elements of 0^ n (u), so that to order 
O(p ) one has 

Bn = f J ^dxx n ~ 1 <p 7T (x) = b n>n -i. (n odd) (53) 

Combining this with l|48|l and going from moments to rr-space, we obtain the soft pion theorem 
H I=1 (x,£ — ► 1,0) = (j>Tr(x) from |14j . Let us now consider the order 0(p 2 ) corrections to the one-pion 
matrix elements. A nonanalytic dependence on the pion mass is generated by one-loop graphs with 
the insertion of the leading-order operator (|50|). which according to (|32j) contains terms with three 
pion fields. Since for different n the operators in Q50|) only differ by their coefficients 6 nin _i and the 
number of overall derivatives ud, the one-loop corrections of the one-pion matrix elements for all n 
are universal and hence the same as the corrections to the isovector axial current, given by (|52|). We 
can therefore write 

(tt 6 (p)| Of (ft (u) |0> = -i5 ab { P u) n F n 6 n , n _! + Oim 2 ) (54) 

to order 0(p 2 ) accuracy, where the 0(m 2 ) terms are generated by tree- level insertions of higher-order 
terms in the ChPT expansion of 0^ n (u), which can be constructed in analogy to those in Section 0J 
As a consequence, the relation (|53[1 does not receive nonanalytic corrections in the pion mass at next- 
to-leading order in ChPT, and one has B n = b n ^ n -\ + 0{m 2 ) up to terms of order 0(p 4 ). Together 



with the 0(p 2 ) corrections of the form factors A^u~(t), this will allow us to discuss the one-loop 
corrections to the soft pion theorem in Section [5J We note that to the order given, the bare pion mass 
m in (|52[) and (|54jl can be replaced with its renormalized value m n , given by |17j 



m 



m 



m 



log + 



2m 2 



(55) 



32n 2 F 2 b fi 2 F 2 

Let us finally show that our representations fUjl and (|4"5|) of the operators O s n {u) and 0^(u) in 
ChPT are equivalent to those given in |19j without derivation. There, the matching was done for the 
nonlocal operators appearing in the definition of GPDs. For the left-handed current we find in our 
notation 

1-75 



-\\u) n^7 M - 

If* 



dp 



l- 



da F(0, a) Uj k (Ai«) (iud ) U kj (A 2 «) + . • • , 



(56) 



-i J-l+\P\ 

where Ai = —\{f3 + a) A and A2 = — a)A, and where the dots denote terms not contributing to 
two-pion matrix elements at order 0(p 2 ). Due to time reversal invariance the function F(/3, a) is even 
in a. Taylor expansion gives the local operators 



1 - 75 



(iuD) n - 



(57) 



F 2 Y ( ) d(3 dap n - m - 1 a m F(p,a) (UudY 



+ 



where all fields are taken at space-time position zero and the dots have the same meaning as before. 
By repeated use of the identity dU = Ud + UL we can rewrite 



(ud) n U 



n-2 

E 

m=0,2,. 
n— 3 

E 

m=0,2,. 



fray 



L(ud 



\n—m—2 



L(ud) 



n—m—2 



L 



L 



+ C(L 3 ) (n even) 



+ {{\ud) n ~ l L + 0(L 3 ) (nodd) 



(58) 
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and then project out the isoscalar part for n even and the isovector part for n odd. Repeating the 
same arguments for the right-handed operators, we recognize precisely the operators that appear in 
(JUJ and (@SJ). 



4 Tensor operators at order 0(p 2 ) 

The number of independent terms in the ChPT expansion of operators rapidly grows with the order. 
Henceforth we shall consider only the terms which contribute to two-pion matrix elements at tree- 
level. Taking into account the discussion in the previous section, one finds that the corrections of 
order 0(p 2 ) to the operators can be cast into the form 

1 1 " /2 

AO s n {u) = - a* tr [x- L n (u) + X - Rn{u)) + - £ a* fc tr [L k {u) L n ^ k {u) + R k {u) R n „ k (u)} tr X + 

k=l 

n/2 

+ -J2^n,k^{L p ^)K-k( u ) + R P A u )K-k( u )} ( 59 ) 
k=0 

for 0^(u) and 

(n-l)/2 

A0» = ibl (L a n (u)+RUu))tix L + + e abc £ 6* fc (l£(u) L« _ k (u) + R b k {u) K_ k (u)) tr X ^ 

fe=l 

(n-l)/2 

+ e abc 2 + (60) 

fc=0 

for C"(u), where we have introduced L P:k (u) = ^p W ...^ fc w w . . . « M * and its right-handed counterpart. 
Using hermiticity and time reversal invariance one can deduce that the coefficients a*, b*, a* k , b* k 
and Aa njk , &b n ^ k are real. The order 0(p 2 ) part of a twist-two operator which contributes to two-pion 
matrix elements is hence parameterized by 2 ([n/2] + 1) real constants. Note that the number of the 
form factors k (t) in Q is [n/2] + 1. At order 0(p 2 ) the form factors thus read 

<*(*) = A n!k + A nf ™ 2 + A n,k * + loo P contributions, (61) 

where all parameters A^'I"^ , A^' k are independent and can be expressed in terms of the coefficients 
a* 6* etc. in © and 



5 One-loop contributions 

To calculate the two-pion matrix element of the operator O s n (u) + AO*(ii) or 0®(u) + AOJJ(«) one 
has to take into account loop contributions. The diagrams to be calculated are drawn in Fig. ^ The 
technique of such calculations in ChPT is well elaborated, so we omit the details and directly give 
our results. Introducing the notation L a (Xu) = u^L^(Au) and its analog for the right-handed fields 
one gets 

m 2 — 9/ 

(vr b (y)|^(Am)L c (A 2 n)K(p)} lloo P = - 6«* IjL_fJ r(e) 



4m 



(l-a 2 )t 



WlTfl 2 



(62) 
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a 



b 



c 



Figure 1: The one-loop graphs contributing to the two-pion matrix elements of the twist-two operators 
O s (u) and O a {u). The black blobs denote the operator insertion and the black boxes the 7r 4 vertex. 



in d = 4 — 2e dimensions, where \x is the renormalization scale. Taylor expanding in Ai and A2 as 
we did in going from 1)56)1 to (|57|). one can read off the one- loop corrections for the local operators in 
(|41j) . Taking into account ©, one concludes that the form factors A^®(t) with k < n do not receive 
additional contributions at one-loop order and hence depend linearly on m? and t, 



a(0) + A ihm) m 2 

A n,k + A n,k m 



A (1,t) t 



(fc < n - 2) 



(63) 



up to corrections of 0(p 4 ), where as in Section we omit isospin labels in 4°L 4'fc and A^'^ . For 
the form factor Af^(t) one finds after some algebra 



i(i.t) 



4=°(i) = 4% + 4 1 ^)( M )m 2 + 4 1 ^( M )* + 



2* 



64vr 2 F 2 



(64) 



^ 2 fe - n A (0) 



fc=2,4,. 



-1 



da (1 — a )a 



n—k 



log 



4m 



(1 - a 2 )t 



4/i 2 



+ 1 



where A\ I is given in 61. Here AX™ fa) and 4^V) 

are renormalized using the subtraction 
scheme of ^7j, where together with a 1/e pole one subtracts a finite constant in the combination 
1/e + V^(2) + log(47r). In particular, we find for n = 2, 4 



,(l,m) 



^2,2 



A 



(0) 



2.2 

(l,m) 



1 + 



m 



2t 



48tt 2 F 2 
• 2 2/ • 



m 2 4 t + 2m 2 . . 

log — + J(t) 

\i z 3 i 



.4 



(0) 
4,4 



1 + 



+ A 



(0) 



40tt 2 F 2 
2 2i 



m" 11 1 9 1 . * , . 



4 - 2 960vr 2 F 2 



1 



log -j- - — + a 1 + —(3a 2 - 5) J(t) 



15 



(65) 



(66) 



where 



J(t) = 2 + crlog 



a 



a + 1 



4m 2 



a 



(67) 



For |t| < 4m 2 one can expand J(t) = -2/(3o- 2 )-2/(5o- 4 )-0(l/cr 6 ) = t/(6m 2 )+t 2 /(60m 4 )-O(t 3 /m 6 ). 
Note that to the order we are working at, one may replace the bare quantities F 2 and m 2 in the 
formulae of the present section with their one- loop expressions, given in 1)52(1 and 1)55)1 . 
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The order 0(p 2 ) correction to the matrix element of the isovector vector current (|2()j) reads JJj 



(ApOl^fKW) =2ze abc P ll (l + F v (t,m 2 )-ll( f ,)t/F 2 ) , 
where the renormalized low-energy constant Zg(/i) is defined in and 



(68) 



F v (t,m 2 



1 



128tt 2 F 2 
1 

96vr 2 F 2 



da 



4m 2 -(l-a 2 )t log 



4m 2 - (1 - a 2 )t 



2 2 

m m 
log- 



l 67T 2 F 2 *> „2 



(t - 4m 2 ) J(i) - i - + log -^r 



(69) 



To complete the calculation of the isovector current at order 0(p 2 ) it is sufficient to compute the one- 
loop corrections to the pion matrix element of the nonlocal operator O a (Xi, A2) = e abc L b {\\u)L c (\2u). 
We obtain 



<7rV)|0 c (Ai,A 2 )KGp)) lloop 
1 



{n b (p')\O c (X u X 2 )\n a (p)) 



a/„\\ tree 



64vr 2 F 4 
1 



16tt 2 F 2 v ; \4ir/i 2 

r(-l + e)£(Pu)V ?(Pu) (Al+Az) 



' / da 



Am 2 - (1 - a 2 ) t 



Am 2 - (1 - a 2 )t 

167T^ 2 



d_ 

da 



3 <af(P«) (Ai-Aa) (1 _ Q 2 )l . (70) 



Making use of (|45|) and (|47j) one can represent the form factors A n ^(t) with k < n — 1 as 



^ (* ) = <i ^ " ^2 log ^ J + (/*) m a + 4ff t (k < n - 3) 

up to corrections of order 0(p A ). For the form factor A^^L-^i) one finds after some algebra 



m 



m 



,(l,m), 



(71) 



n-1 



(72) 



+ 



1 



64^ E 2 fc -(n-^4°L 



fc=0,2 

For n = 1 this gives 
i(i,t) 



da a 



n— fc— 1 



4m 2 -(l-a 2 )i log 



4m 2 - (1 - a 2 )t 



A^t) = 2(l + F v (t, m 2 )) + A^ } (p) t , 



(73) 



where A\q (/x) = —21q(p)/F 2 according to l|68|). Note that at i = this form factor is fixed by 



current conservation, so that A^Q L \pb) = 0. For n = 3 we find 

4f(t) = Afl(l + F v {t,m 2 ))+A { if^)m 2 + A^^)t 



A 



(0) 

3,0 



t 



1280vr 2 F 2 



9 v , ?n 2 6 4 9 , . 

(5a 2 - 3) log — + - - -a 2 - 2^J{t) 



(74) 



2 The corresponding form factor is known to two-loop accuracy now |2()j . 
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Our results (|7T|) and (|72*|) agree with Eq. (27) in ^H] , and our results and (|6*4"|) agree with Eq. (26) 
there if a misprint is corrected. 3 Furthermore, our expression (|65j) for -A^ ^) a g rees with |16| . 

From ()63j) and ()64j) we can readily generalize to order 0(p 2 ) the relation (|44|) for the sum of 
isoscalar form factors of given n. In particular we find for t = 

n 2 2 ^ qj> 

^ E ^(0) = jj&j log ^ ^ E_ ^(0) („ + 1)( f_ t+1) + °«>> PS) 

where on the right-hand side we have replaced the tree-level quantities F, m, A^ k by the full ones. 
We have a relation between observables, which may for instance be tested in lattice calculations as a 
consistency test for the applicability of one-loop ChPT for a given range of unphysical pion masses. 
Using Q71|) and (|72|) and our discussion following (|54jl. we obtain a relation at order 0(p 2 ) between 
the sum of isovector form factors of given n and the moments B n of the pion distribution amplitude. 
For t = the logarithms in m 2 cancel and we have 

n-l 

£ 2 k A I = k 1 (0) = rB n + O(ml). (76) 

fe=0,2,... 

For n = 1 this relation is trivial since J 4{^ 1 (0) = 2 and B\ = 1 are exact, but for re > 3 it presents 
again an interesting test in lattice QCD, where both the form factors A^ k (t) and the moments B n 
are calculable. Translated to x-space, one finds that the soft-pion theorem H I=1 (x,£, — > 1,0) = 
<fin(%) + C( m -w) does not contain nonanalytic terms in the pion mass at one-loop order, as was pointed 
out in [in]. 4 



6 Summary and conclusions 

We have calculated the moments of the generalized parton distributions of the pion at order 0(p 2 ) of 
ChPT. Our main results are given in ()63|) . (|64|) and (|71jl. (|72|) for the isoscalar and isovector matrix 
elements, respectively. We find that at one-loop order only the form factors accompanying the highest 
power of £ at given n have a nontrivial dependence on t and m 2 . All others depend linearly on t. The 
isoscalar form factors are linear functions of m 2 as well, whereas the isovector ones receive universal 
logarithmic corrections in m 2 . At order 0(p 2 ) the [n/2] + l form factors A l n k at given n are determined 
by 3 ([n./2] + 1) real parameters, which in the isoscalar case are subject to the restriction (|44|). One 
can check that for a fixed ratio t/m 2 the form factors in 1)64(1 and (|72[) deviate from a linear behavior 
in only for pion masses close to the physical value, and that this deviation is rather small. 

The possibility to use the formulae derived here to extrapolate lattice data obtained for pion 
masses around 500 MeV to the physical pion mass obviously depends crucially on the size of the 
higher-order corrections. Taking into account that ChPT works quite well in the kaon sector [22] one 
may hope that the two-loop corrections will not be large for m 2 ~ 500 MeV. In any case a strong 
deviation of the lattice data from a linear dependence on m 2 would indicate the need for higher-order 
ChPT fits to describe the lattice data. The comparison of the ChPT predictions with lattice data 
for the form factor A^ 1 ^) = 2F£ m (t), where all 0{p 2 ) parameters are known |17M20| . will provide a 
crucial convergence test. The same holds for relations which involve only lattice observables, namely 
the relation (|75|) between isoscalar form factors and the relation l)76|) between isovector form factors 
and the moments of the pion distribution amplitude, as well as their generalizations for finite t. 

3 The one-loop term in Eq. (26) of needs to be multiplied with a factor 1/2 21 . 

4 To be precise, in |19| this relation was given with the pion distribution amplitude taken in the chiral limit, without 
discussing its one-loop properties in ChPT. 
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At two-loop order six new parameters appear for each form factor. Three of them are the co- 
efficients of the analytic terms tm 2 , t 2 , and the other three arise from terms in the expansion 
of operators which do not contribute to the two-pion matrix element at one- loop order (m^ log 2 m\ , 
m^logm 2 , and tm\ logm 2 ). Obviously the lattice data, even if it covers several different volumes 
and many pion masses, can hardly be hoped to be accurate enough to obtain a stable simultaneous fit 
of all these parameters. However, the same parameters should enter the m T and volume dependence 
of other lattice observables, such that these will provide complementary information. A detailed 
analysis of these relationships remains, however, to be worked out. 
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